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I. INTRODUCTION 



The process of photo-recombination in highly charged heavy ions has been the subject of many 
theoretical and experimental studies up to today, concerning both radiative recombination (RR) 
and dielectronic recombination (DR) (see, e.g., tlL|2j]) and their interference. With the enhanced 
experimental possibilities and achieved precision, the subject of electron recombination into highly 
charged ions has been expanding to include QED corrections yfl. The effect of interference be- 
tween RR with DR has been theoretically studied (see, e.g., Ref. J4]) and experimentally con- 



cluded D2L 



Q. 



In Ref. [|7|] a recombination process that is the nuclear analog of DR has been theoretically pro- 
posed. Although not yet experimentally observed, nuclear excitation by electron capture (NEEC) 
has been an interesting subject after experimental observations of atomic physics processes with 
regard to the structure of the nucleus have been recently reported, such as bound-state internal 
conversion [8|] and nuclear excitation by electron transition (NEET) [9]. In the resonant process 
of NEEC, a free electron is recombined into a bound state of an ion with the simultaneous exci- 
tation of the nucleus. The excited nucleus can then decay radiatively or by internal conversion. 
Several theoretical studies have been made concerning NEEC in plasmas [3, [k]] or in solid tar- 



gets 







1311 . In 111411 we presented relativistically correct theoretical cross sections for NEEC 



followed by the radiative decay of the nuclear excited states for highly charged heavy ions. 

If the initial and final states for NEEC and RR coincide, quantum interference between the 
two processes occurs. Such an interference effect is interesting as it involves two very different 
pathways: while in RR only the recombining electron is involved, NEEC corresponds to a quan- 
tum path in which the nucleus is excited. In Figured] the RR and NEEC mechanisms are shown 
schematically. Beside NEEC, the strong competing process of RR is always present in an exper- 
iment. Therefore, the magnitude of the interference effect may also play an important role for 
observing NEEC. 

In this paper we theoretically investigate the interference between NEEC and RR, focusing 
on collision systems with suitable excitation energies that could be candidates for experimen- 
tal observation. We derive the total cross section of the recombination process with the help of a 
Feshbach projection operator formalism, which allows the separation of the interference term from 
the NEEC and RR cross sections. The radiation field is expanded in terms of multipoles in order 
to clearly discern the NEEC transition multipolarities. The electric and magnetic electron-nucleus 



2 




FIG. 1: NEEC and RR recombination mechanisms of a continuum electron into the K shell of a bare 
ion. The nucleus is schematically represented as undergoing the transition from the ground state (G) to the 
excited state (E) and again to its ground state. 



interactions are considered explicitly, and the nucleus is described with the help of a geometrical 
collective model and making use of experimental data. The dynamics of the electron is governed 
by the Dirac equation. We express the interference term of the cross section using the dimension- 
less Fano profile parameter for electric and magnetic transitions in Sec. II. The numerical results 
of the calculation are given in Sec. Ill, together with an interpretation of the results regarding the 
possibility of an experimental observation of NEEC. We conclude with a short summary. In this 
work atomic units are used unless otherwise specified. 



II. THEORETICAL FORMALISM FOR INTERFERENCE EFFECTS 

In this section we derive the total cross section of the recombination process involving NEEC 
followed by the radiative decay of the nucleus and RR by means of a Feshbach projection operator 
formalism. We consider that the electron is captured into the bound state of a bare ion or an ion 
with a closed-shell configuration. We calculate the interference term between NEEC followed 
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by the radiative decay of the excited nucleus and RR in the total cross section for electric and 
magnetic multipole transitions of the nucleus. 

A. The interference between RR and NEEC in the total cross section 

The initial state j^) of the system describing the nucleus in its ground state, the free electron, 
and the vacuum of the electromagnetic field can be written as a direct product of the nuclear, 
electronic, and photonic state vectors: 

|*i> = \NIiM Iv pm s ,U) = \NUM h ) ® \pm s ) ® |0> . (1) 

Here, p is the asymptotic momentum of the electron, m s its spin projection, and \N) the nuclear 
ground state, denoted by the total angular momentum Jj and its projection Mj. . In considering RR 
or NEEC followed by the radiative decay of the nucleus, the final state of the recombined system 
\tyf) consists of the nucleus in its ground state, the bound electron and the emitted photon. Rather 
than using the plane wave expansion for the electromagnetic field as in J14I1 . it is more convenient 
in this case to consider photons of a given angular momentum and parity. The final state can be 
written as 

|*/) = \NI f M If ,n f K f m f ,XkLM) (2) 
= \NI f M If ) <g> \n f K f m f ) <g> \XkLM) , 

with rif, Kf, and mj being the principal, Dirac angular momentum, and magnetic quantum num- 
bers of the bound one-electron state, respectively. The emitted photon has the wave number k, 
the total angular momentum L and its projection M. Furthermore, A stands for electric (e) or 
magnetic (m) waves. The intermediate resonant state formed by the electron capture in the pro- 
cess of NEEC consists of the excited nucleus, the bound electron, and the vacuum state of the 
electromagnetic field, 

\V d ) = \N*I d M Id ,n d K d m d ,0) (3) 
= \N*I d Mj d ) <g> \n d K d m d ) g> |0) . 

The excited nuclear state is denoted by \N*). In our case, the recombined electron does not 
undergo further decay cascades, i.e., n d = nj, K d = Kf and m d = nif. 
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Following the formalism presented in Ref. 111411 . we introduce projector operators onto the in 



dividual subspaces, in order to separate these states in the perturbative expansion of the transition 



operator. We neglect corrections due to two- or more-photon states H15L I16I1 and due to the presence 
of the negative electronic continuum. The Fock space is then given by the sum of three subspaces: 
the subspace of the states that contain the free electron, with its projector operator P, the subspace 
of the states characterized by the presence of the excited nucleus, together with the correspond- 
ing projector operator Q, and finally the subspace of the states with a photon, associated to the 
projector operator R. We postulate the completeness relation 

P + Q + R = 1 , (4) 

where 1 is the unity operator of the total Fock space. 

The total Hamiltonian operator for the system consisting of the nucleus (n), the electron (e), 
and the radiation field (r) can be written as 

H = H n + H e + H r + H en + H er + H nr . (5) 



The expressions of the first three Hamiltonians can be found in UM- Interactions between the three 
subsystems are described by the three remaining Hamiltonians in Eq. ©. We adopt the Coulomb 
gauge for the electron-nucleus interaction (en) because it allows the separation of the dominant 
Coulomb attraction between the electronic and the nuclear degrees of freedom: 

H en = [ d\ n p^- . (6) 
J \r e -r n \ 

Here, p n (r n ) is the nuclear charge density and the integration is performed over the whole nuclear 
volume. The static part of the electron-nucleus interaction is contained in Hamiltonian H e . The 
interaction of the electron with the transverse photon field quantized in the volume of a sphere of 
radius R is given by 

H er = -a ■ A = - ^2 ( a lfcLAf« ' A\kLM(r) + H.c.J , (7) 

XkLM 



with the vector potential of the quantized electromagnetic field Ill7h 



Ar)= (^kLM{r) a\ kLM + A* XkLM (r) a XkLM ^j ■ (8) 



XkLM 
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Here, a is the vector of the Dirac matrices and the two independent solutions of the wave equation 
for the AxkLMir) are 



XnfiLMlF) = \l^h(kr)Y*l(0,<p), (9) 



where the quantum number k is discretized by requiring the proper boundary conditions at a per 
fectlyc 
by 111 



fectly conducting sphere of radius R. The YH(6, cp) denote the vector spherical harmonics, given 



(e x + ie y ), (11) 



Y%(e, (p) = E C ( L 1 L; v v M ) Y ^o, v% , do) 

V q 

where q — 0, ±1 and the spherical unit vectors e q expressed in terms of the Cartesian unit vectors 

i^xi ^yi &z) 

1 

e*o = e 2 , 

e_ = —j=(e x — ie y ) . 

Similarly, the interaction of the nucleus with the electromagnetic field is given by the Hamilto- 
nian 

H nr = — ^2 \ a \kLM I d 3 r n j n (r n ) ■ A Xk LM(r n ) + H.c.j , (12) 

C XkLM \ J / 

where j n (r n ) is the nuclear current. 

Using the projection operators we can separate the perturbation V in the total Hamiltonian 

H = H + V, (13) 

with 

H = PHP + QHQ + RHR , (14) 
V = H-H = PHQ + QHP + PHR (15) 
+ RHP + RHQ + QHR. 

This way of defining H has the advantage that the effect of the nuclear potential on bound and 
continuum electron states is included in H to all orders. The individual terms in the perturba- 
tion operator describe transitions between the different subspaces. For example, QHP describes 



in the lowest order the time-reversed process of internal conversion (IC), namely, NEEC, while 
PHR and RHP are the first-order operators for photo-ionization and radiative recombination, 
respectively. 

n 

The transition operator is defined as LI 9] 

T{z) = V + VG(z)V , (16) 
where G(z) is the Green operator of the total system given by 

G{z) = (z - H)- 1 . (17) 

Here, z is a complex energy variable. The total cross section for a process can be expressed by 
the modulus square of the matrix element of the transition operator, after summing over the final 
states and averaging over the initial states that are not resolved in the experiment, 

»w<*) = f E E 2p77+T) E 08) 

M If m d XLM y 1 'M^m,, 

i- J dn p \hn + \{* f \T(E + ie)\%)\ 2 Pf , 

with the \& f and ^ as final and initial eigenstates of H , respectively [see Eqs. ((U) and ©]. Here, 
Fi denotes the flux of the incoming electrons, pj the density of the final photonic states, and VL p is 
the direction of the incoming free electron characterized by the angles 9 P and <p p . 
We use the Lippmann-Schwinger equation 

G(z) = G {z) + G (z)VG (z) (19) 
+ G (z)VG (z)VG (z) + ... 

to write the perturbation series for T(z) in powers of V with the Green function Gq(z) of the 
unperturbed Hamiltonian H : 

T(z) = V + VG (z)V + VG (z)VG (z)V + ... . (20) 

Since the initial state of the NEEC process is by definition an eigenstate of P, and the final state is 
an eigenstate of R, we only need to consider the projection RTP of the transition operator: 

RT(z)P = RVP + RVG (z)VP (21) 
+ RVG (z)VG (z)VP 
+ RVGo(z)VG (z)VG (z)VP + . . . . 
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The first term RVP accounts for the radiative recombination. Taking into account from the infinite 
perturbation expansion in Eq. (f2TT) the terms that correspond to NEEC Ill4ll we can write the final 
expression for the transition amplitude for the recombination process as 

(V f \RT(z)P\%) = (VflRHerPlVi) (22) 

^ z-E d + ^Y d 

Here, Y d denotes the total natural width of the excited state \d) = \N*I d M Id ,n d K d m d ,0). The 
magnetic interaction Hamiltonian H magn accounts for the recombination of the free electron by 
exchanging a virtual transverse photon with the nucleus in the unretarded approximation [14], 

H m ann & [ d T n jz^ ~ 7 ■ (23) 

c J \r-r n \ 

Using the expression of the transition operator, the total cross section can then be written as 

^ - 1 E E ^ Es/«- ^ 

1 M If m d XLM y % > Morris J 

| (NI f M If , n d K d m d , XkLM\H er \ NliMj. , pm s , 0) 



(E-E,) + lT d 
x (N*I d M Id , n d K d m d , 0\H en + H magn \NIiM h ,pm s , 0) 



Pf 



The first term in the modulus squared accounts for RR and the second one for NEEC. We can 
separate therefore the equation above in three terms, 

a^ f {E) = a RR (E) + a mEC {E) + a- mt (E) , (25) 

with the RR and NEEC total cross sections given by 

ffRR(E) = | £ £ ^ E (26) 

1 M If m d \LM y 1 ' Mj.m 3 

J 1 

/ dQpliNIfMT^ndKdm^XkLMlHerlNIiM^pms^^pf, 



and 



1 M If m d M Id \LM y 1 ' M h m a 



— I (lil, 

47T 



(NIfMi ,n d K d m d , XkLM\H nr \N*I d M Id , n d K d m d , 0) 



(E - E d ) + f T d 



x (N*I d M Id ,n d K d m d , 0\H en + H magn \NIiM h ,pm s , 0) 



PS 
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The term describing the interference between RR and NEEC can be written as 

1 M If m d M Id XLM V 1 ' Morris 

J_ /" / (NI f M If ,n d K d m d , XkLM\H nr \N* I d M Id , n d K d m d , 0) 

4nJ p { (E — E d ) + |r d 

x (N*I d M Id , n d K d m d , 0\H en + H magn \NIiM h ,pm s , 0) 

xiNIfMj^ndKdm^XkLMlHerlNIiM^pm^O)* + H.c.J . 

The aim of this paper is to calculate the interference term in the total cross section. The calcu- 
lation of the NEEC cross section and predicted values for several collisions systems can be found 



in jl4\ . Furthermore, the calculation of the RR total cross section is well understood. An extensive 



tabulation of relativistic total cross sections for RR as a function of energy ranging from closely 



above the threshold to the relativistic regime of relative electron energies is available in |20(]. 

If we consider the matrix element of the Hamiltonian H er connecting the radiation field and the 
electrons in the interference term, the initial and the final nuclear total angular momenta as well as 
their projections have to coincide, as they are not influenced by RR, 

(NI f M If ,n i K d md,MLM\H er \NI i M Ji ,pm.,0) (29) 
= Si^.Smj.Mj. (n d n d m d , XkLM\H er \pm s , 0) . 

The initial state continuum electronic wave function is given through the partial wave expan- 



sion 021D 



\pm s ) = ^ AK Y; mi (Q p )C (l X - 3 - mi m s m) \sKm) , (30) 

Kmmi ^ * 

where e is the energy of the continuum electron measured from the ionization threshold, e = 



^/p 2 c 2 + c 4 — c 2 . The orbital angular momentum of the partial wave is denoted by / and the 
corresponding magnetic quantum number by mi, while the partial wave phases A K are chosen so 
that the continuum wave function fulfills the boundary conditions of an incoming plane wave and 
an outgoing spherical wave. The total angular momentum quantum number of the partial wave is 
j = | k | — ~ . The interference cross section in the case of NEEC involving a nuclear transition with 
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specific parity A and multipolarity L can then be written as 



'hit 



*l v v ei y- (3d 

Ft ^ ^ 2(2 Ji + 1) ^ 4vr 



M Id M T . Mm d 



( (NI.M^ XkLM\H nr \N*I d M Id , 0) „ . 

x s-T- — j- (JV JdM/j , n d K d m d | H en + H magn \ NIiM h , sum) 

V \E - E d ) + 2 1 d 

x (n d K d m d , \kLM\H er \sKm, 0)* + H.c. 

We can relate the interference cross section term with the NEEC cross section, introducing the di- 
mensionless Fano profile parameter Qf. The expression of the NEEC cross section from Ref. Ill4ll 
is 

2tt 2 Ad-*fY i -^ d 

a^ f (E) = — r n L d (E - E d ) , (32) 

v r d 

where Af^f is the radiative rate defined as 

2tt 



At"** 



and is the NEEC rate, 



- Yl Y\{NI } M Ir n d K d m d ,\kLM\H nr \N*I d M li ,n d K d m d M 2 pf (33) 



^ d ' " .W, M Mi 

f a 



Yt d = 2f2/ 27T n Yl Yl dn p \{N*I d M Id ,n d K d m d ^H en +H^^ 

^ 1 ' M h m s M Id m d J 

(34) 

Furthermore, p denotes the continuum electron momentum and p, t the density of the initial elec- 
tronic states. The explicit energy dependence of the interference term can be expressed with the 
help of the Lorentz profile L d (E — E d ), defined as 

Ld{E - E d ) = ^Z^TFI » < 35 > 

which in turn is related to the NEEC total cross section. The interference cross section can be 
written in the concise form |4J] 

= CTneeo i^1ttt ( 2£ ^r Re H) + Im (i) ) • <36) 

with the dimensionless Fano profile parameter 

= Y J~] y~](N*IdMi A , n d K d m d \H en + H magn \NIj M Iv ektu) (37) 

MijMi. Mm d Km 

(N. 

x 



(NliMii, XkLM\H nr 


N*I d M Idl 0){n d K d m d , XkLM\H er \t 


ZKm, 0)* 


H I d 


(NliM'j., kXLM' 




N*I d M' d ,0)\ 2 
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We have used prime indices for the summations in the expression of the nuclear radiative rate in 
the denominator in order to avoid confusion. With the further observation that the Fano profile 
parameter l/Qf is real for both the electric and magnetic cases, the interference cross section 
yields 

2(E - E d ) 2I d + 1 1 

*i=t - ^EEC ^ 2Ii + 1 Q- f ■ OS) 

B. Electric transitions 

In order to calculate the matrix elements in the Fano profile parameter in Eq. (|37l) , an adequate 



nuclear model is needed. Following the outline in II 1411 . we describe the nucleus by means of a 



geometrical collective model [22] which assumes that the excitations of the nucleus are vibrations 



and rotations of the nuclear surface, which is parameterized as 

oo t 

R(9, <p, t) = R (l + E a im(t)Yim{B, <p)) . (39) 

1=0 m=-l 

The time-dependent deformation amplitudes at m (t) describe the nuclear surface with respect to a 
sphere of radius R and serve as collective coordinates. This parameterization can be used to cal- 
culate the matrix element corresponding to the NEEC process for a given partial wave component 



and a given multipolarity L, that yields 111411 



{N* hM^n^m^H^NUM^eKm) (40) 

L 

^ ( _l ) / d +M, I+ L +M+m+ 3 Jd ^-(L + 2)^^^ K(iv , /(i||gL||iv/j) 



fi=—L 



x y/2j d +lJ (2L ^ 1)3 C(^ h L] -M u M Id ^)C{j j d L; -m m d - fiC (j d Lj;^0 ~ J 



where Q LM is the electric multipole moment defined by Ill7ll 



Qlm = J d 3 r n r^Y LM (9 n ,Lp n )p n (r n ) . (41) 



The electronic radial integral is given by 

-Ro 



RL, KdlK = / drrL+2 (fn dKd (r)fe K (r) + 9n dKd (r)g £K (r)) + (42) 



roc 

- < +2 / drr- L+1 (f ndKd (r)f £K (r) + 9n dKd {r)g £K {r)) . 

'Ro 
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with g £K (r) and f eK (r) being the large and small radial components of the relativistic continuum 
electron partial wave function 

.-.(0=f «~ WI * W "* J V (43, 

with the spherical spinor functions fi™, and g nd K d (r) and f nd K d (r) the radial components of the 
bound Dirac wave function. 

For the matrix element of the interaction Hamiltonian (PT21) between the nucleus and the radia- 



tion field, we follow the outline in |23j], considering that the wavelength of the radiation is large 



compared to the nuclear radius, kR <C 1, so that the Bessel functions can be approximated in the 
first order in kr as 

^ =< ' <44) 

In this case the electric solution of the wave equation can be written as 



A [e)kLM (r) = -y — ^- - — (2L + 1)\\ ^ ' Y LL-A e i <P) ■ (45) 

With the use of the continuity equation for the nuclear current j n we obtain for the matrix element 



(NIiM h ,(e)kLM\H nr \N*I d M Id ,0) = (-1)^-^+1 \j —C (I, h L; M h -M h M) (46) 



y/L(2L + l) (2L + 1)!! V '" VL " d/ 
The remaining matrix element of H er can be evaluated by writing the electric solution of the wave 
equation in Eq. © in a more suitable form. Using the properties of the vector spherical harmonics 



Il24ll we obtain 



4:7Tck I I L , "Af 



A (e)kLM (?) = \^^—l3L + i{kr)Yt L+l {e^) (47) 



L + 1 "iwWS- 1 (^) 



2L + r 

The electron-radiation interaction matrix element then yields 

(n d K d m d , (e)kLM\H er \eKin, 0) (48) 



'"^ ^ (n d K d m d \j L+ i(kr)a ■Y^[ +1 (9,Lp)\eKm) 



R V V 2L + 1 



'n d K d m d \j L ^ x {kr)a ■ Yjf L _ x {0 , ip)\eKm) 



2L + 1 
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The matrix elements containing the product of the Bessel spherical functions, the Dirac matrix a 
and the vector spherical harmonics can be expressed in a compact way using the properties of the 
spherical tensor operators Il25ll . The expression in the above equation becomes 



(ridKdfnd, (e)kLM\H er \eKm, 0) 



(49) 



i(-ly + M ——C{j L j d ;m M m d 



R 



4tt 



2 2 u 



X 



L + l 



+ 



L(2L+T) (LJ£Ll " (Kd "' 6)itl) 

<(l + 1) + (« d - k) 



(L + l)(2/ + l) 



with the radial integrals 



drr 2 j L (kr) (g niKd (r)f £K (r)±g £K (r)f ndKd (r)) . 



(50) 



Combining the formulas of the three matrix elements from Eqs. (l40l) . (|46|) and (|49|) in the expres- 
sion of the Fano profile parameter Qj and using the summation properties of the Clebsch-Gordan 
coefficients we obtain the final formula 



1 



(e) 



i?-^(2 Jd + i; 



(L + 1)(2L + 1) 



(51) 



fc- £ (2L+l)!!j3fl w (2j + l) ( ^ \ L 



2 



L + l 



+ i 



1 -((L + l)/ L+1 + (^- K )/+ +1 ) 



(L + 1)(2L + 1) 



C. Magnetic transitions 



The magnetic transitions in the nucleus can be easily included in the calculation by assuming 
that the electron does not penetrate the nucleus, i.e., that the electronic radial coordinate r e > r n 
is always larger than the nuclear radial coordinate. This approximation is precise enough for the 



studied cases [2 



3, 



27|]. The NEEC matrix element for the magnetic transition, involving only the 
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magnetic Hamiltonian H magn for a given partial wave and a given multipolarity can be written 



as ||14il 



{N*I d M Id ,n d K d m d \H magn \NIiM h ,eKm) = (52) 
4 "\/ l(2L + + 1 1)3 X>1) 7 '~ M/ ' +m+1 C(I d h L; M Id - M h ii){-N*I d \\M L \\NI % ) 
x (n d K d m d \r~ {L+1) a ■ Y^(9, ip)\eKm) , 



where the electronic matrix element can be evaluated in a similar way as the ones in Eq. (1491 ) to 
yield 

(n d K d m d \r-^a-Y^(e,(p)\eKm) = (53) 



oo 



x / drr L+1 (g ndKd ( r )fe.(r) + f nAKd (r)g £K (r)) . 



1 2 





This way of writing the electronic matrix element is equivalent to the more lengthy one presented 
previously in [| 1-4T | . 

Now let us consider the matrix element corresponding to RR. It has, up to the presence of the 
spherical Bessel functions, a similar expression, 



(n d K d m d , (m)kLM\H er \eKm,0) = -\j —^-(n d K d m d \j L (kr)a ■ Y^ L {e,ip)\eKm) . (54) 



Using the properties of the spherical tensor operators Il25h . we can write the RR matrix element as 



(n d K, d m d , (m)kLM\H er \eK,m,0) = (55) 



i(-l) J L i\ — —-C(iL] d ;mMm d )(K d + K) 

R v ) y 4ttL(L + 1) U Jd ' dA d ' 

X I / drji(kr) (g nd K d (r)fen(r) + fn d K d (r)g eK (r)) ■ 



The remaining matrix element involved in the expression of the Fano profile parameter Qf is that 
of the interaction between the nucleus and the radiation field (fT2)) . We make use again of the 
long-wavelength approximation, so that the spherical Bessel functions are written as in Eq. (144]) . 
With this approximation and using the properties of the vector spherical harmonics, the magnetic 
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solution of the wave equation can be expressed as 



A [m)k LM(^ = \l „ , 7 == = — T , lN „ (r x V)(r L Y LM (6,<p)) . (56) 



' Anck k L 1 
R i^L(L + T){2L + l)\\ 

Rewriting the Hamiltonian H nr we obtain 



Anck L + 1 k If - 

Hnr = l^T\l L (2L + l)!!c(L + l) J dir ^ x 3n{r n )) -V{r^Y LM {e n ^ n )) . (57) 

The integral over the nuclear coordinate can be related to the magnetic multipole operator M LM , 



defined as Ill7ll 



M LM = lT , ,x / rf a r n (f n x j n (r n )) ■ V(r^Y LM {9 n , <p n )) . (58) 



c(L + l] 

The matrix element of the interaction Hamiltonian between the radiation field and the nucleus 
yields 

{NhM h , (m)kLM\H nr \N*I d M Id , 0) (59) 

= h I^^^E±L/ NI . Ml \M LM \N*I d Mn) 
V R v^L(2L + l)!! v 1 /J LMl d Idl 



: _ i)Id _ Midi j4nck k L I L + l 



R (2L + 1)!! y L(2L + 1) 
xC(I d I t L;M Id -M h -M)(NI i \\M L \\N*I d ) 



Combining the results from Eqs. (|52|) . (1551) and (l59l) we write the expression of the dimensionless 



Fano profile parameter Q f , making use of the summation properties of the Clebsch-Gordan 



coefficients: 

1 _^(-D^ +1 (^ + l) r , (2L + 1)!! (60) 



Q f) L(2L+1)(L+1) 

x Y^{2 3 + l)(n d K d + Kf J~ ^rfrr- L+1 ( fed (r)/ £K (r) + / ntiKd (r) fe (r)) ( ^ ^ 

/■oo 

x / drj L (kr) (g ndKd (r)f £K (r) + /„ dKd (r)# £K (r)) . 



/ 

2 

2 -I o 







III. NUMERICAL RESULTS 



We have calculated the Fano profile parameter and the interference cross section term <7 int as 
a function of the incoming electron energy for several collisions systems involving electric E2 
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and magnetic Ml transitions. We consider suitable cases of isotopes which have energetically 
low-lying nuclear levels which make the interference between NEEC and RR possible. 

For the case of the electric transitions we consider the + — > 2 + E2 transitions of the g^U, 
If U, |f Cm, ^ 4 Yb, J|°Er, if Gd, If Gd, If By and If By even-even nuclei. The energies of the 
excited nuclear levels E exc as well as the reduced transition probabilities B(E2), that are needed 
for the calculation of the natural width of the nuclear excited state and the NEEC cross section and 



rate, are taken from Ref. [12811 . The natural width of the nuclear excited state is considered to be 



the sum of the partial radiative rates Af~*f and the IC rates Af c , 

/ i 

Here we sum the radiative transition rates to all possible final states (note that in our case there is 
only one nuclear final state, namely, the ground state). By summing over i we account for internal 
conversion to the initial state of the NEEC process and all other possible IC channels, for the case 
when the capture occurs into a He-like ion. The IC rate can be related to the NEEC rate through 
the principle of detailed balance, 

A d-H _ 2(27, + 1) d 

The NEEC rates and cross sections are calculated using an improved version of the computer 



routines applied in 1141) . We consider the capture into the bare ions of gg Dy, g§ Er, ^ Yb and 



g| 4 Gd. For the cases of the U isotopes and for gg 8 Cm, the capture into the K shell is not possible 
due to the low energy level of the first excited nuclear state. For these three systems, recombination 
into the L shell of initially He-like ions is the most probable one. We regard the capture of the 
electron into a closed shell configuration as a one-electron problem, without the participation of 
the K- shell electrons. We also consider the capture of the electron into the He-like ions of gipGd 
and gg 2 Dy, in which case the width of the nuclear excited state in Eq. (I6TT) contains partial IC rates 
accounting for the possible IC of the i^-shell electrons. 

A numerical evaluation of the radial integrals corresponding to NEEC [RL, Kd ,K, see Eq. (|42l) l 
and the ones corresponding to RR Ul±v Eq- (l50l) l is needed for the calculation of the Fano profile 
parameters and for the interference cross sections. We consider Coulomb-Dirac wave functions 

o 

for the continuum electron and wave functions calculated with the GRASP92 package [|29D by 
considering a homogeneously charged nucleus for the bound electron. In the case of recombination 
into the He-like ions we assume a total screening of the nuclear charge for the continuum electron, 
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i.e., we use Coulomb-Dirac functions with an effective nuclear charge Z e s = Z — 2. For the 
bound electron wave functions, the electron-electron interaction is accounted for in the Dirac-Fock 
approximation. The value of Rl iK , k is not affected by finite nuclear size effects on the accuracy 
level of our calculations. Nevertheless, the finite size of the nucleus has a sensitive effect on the 
energy levels of the bound electron. The energy of the bound electronic state is calculated with 
GRASP92 and includes one-loop one-electron quantum electrodynamic terms, and in the case of 
many-electron bound states approximate QED screening corrections. The nuclear radius Rq is 



calculated according to the semi-empirical formula fl3Qt1 



R = (1.0793A 1/3 + 0.73587) fm , (63) 

where A is the atomic mass number. Values of the Fano profile parameters, as well as the NEEC 
rate and natural width of the nuclear excited state are presented in Table H The values of the 



resonance strength of NEEC, given in ||14fl. 



ad = —r, ^ , (64) 



V r 



(I 



are also presented. 

The Fano line profile parameter characterizes the strength of the interference effects between 
the two recombination channels. Smaller values of \Qf\ indicate more pronounced interference. A 



more quantitative measure of the interference is defined in Ref. Bill as the ratio of the interference 



term and the resonant process term at the energy e±i/ 2 = E d ± T d /2, 



R 



int 



f int (^il/2) 



cneec(£±i/2) 



r ' 2L + l 1 (65) 



rr* 2u + 1 \q,\ ' 



Values for this line asymmetry parameter R mt are given in the last column of Tables HI and IHI1 

A possibility to cross-check the numerical accuracy of the present calculations is given by the 
matrix element of the interaction Hamiltonian H er , which enters the expression of the Fano profile 
parameter. We can use the matrix element to calculate the total cross section for RR for a given 
energy, which can be written in the spherical wave approach as 

c r RR=^-^5^^ / d ^p/^2 \( n d K dni d , XkLM\H er \pm s ,0)\ 2 p f . (66) 

1 m s " rrirf XLM 

RR cross sections calculated by this formula and with the radial wave functions described above 
reproduce the values tabulates in Ref. Il20n with a typical relative accuracy of about one per thou- 
sand, as it can be seen in Table HB 
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TABLE I: Parameters of the NEEC total cross section and the interference term for various heavy ion 
collision systems involving electric quadrupole transitions. E cxc denotes the nuclear excitation energy, E c 
is the continuum electron energy at resonance, Y n stands for the resonant recombination rate, and is the 
total width of the excited nuclear state. The column denoted by S contains the NEEC resonance strengths, 
1/Qf is the inverse Fano line profile parameter, and R mt stands for the profile asymmetry parameter. See 
the text for further explanations. 



Isotope 


^exc(keV) 


£ c (keV) 


Orbital 


Y n {l/s) 


r rf (eV) 




5(b eV) 




1/Qf 


R int 


6 l 4 Dy 


73.392 


10.318 


lSl/2 


1.86 


X 


10 8 


4.37 


X 


10~ 


-8 


3.88 


X 


10" 


-2 


-2.11x10" 


_3 


3.67x10" 


_3 


68 7 °Er 


78.591 


11.350 


lSl/2 


2.23 


X 


10 8 


5.75 


X 


10~ 


-8 


4.70 


X 


10" 


-2 


-2.07x10" 


_3 


4.05x10" 


_3 


ft<Yb 


76.471 


4.897 


lSl/2 


1.79 


X 


10 8 


4.85 


X 


10~ 


-8 


9.27 


X 


10" 


-2 


-2.09x10" 


_3 


4.30x10" 


_3 




123.071 


64.005 


lSl/2 


5.69 


X 


10 s 


2.51 


X 


10~ 


-7 


2.91 


X 


10" 


-2 


-2.61x10" 


-4 


8.77x10" 


-4 


6 fGd 


88.966 


74.742 


2^1/2 


3.35 


X 


10' 


1.21 


X 


10" 


-7 


7.09 


X 


10" 


-4 


-6.10x10" 


_5 


1.67x10" 


_3 


efGd 


88.966 


74.896 


2Pl/2 


1.16 


X 


10 8 


1.32 


X 


10" 


-7 


2.25 


X 


10" 


-3 


-1.16x10" 


_5 


1.00x10" 


_4 


ifGd 


88.966 


75.680 


2p 3 /2 


1.59 


X 


10 8 


1.27 


X 


10" 


-7 


3.17 


X 


10" 


-3 


3.06x10" 


_ 1 


1.86x10" 


Q 
'> 


l?Vy 


80.660 


65.432 


2^1/2 


2.81 


X 


10' 


9.39 


X 


10" 


-8 


6.25 


X 


10" 


-4 


-1.28x10" 


_4 


3.26x10" 


_3 


l?Vy 


80.660 


66.594 


2pi/ 2 


1.59 


X 


10 8 


1.11 


X 


10" 


-7 


2.98 


X 


10" 


-3 


-5.78x10" 


-5 


3.06x10 


-4 


162r) v 

66 u y 


80.660 


66.492 


2P3/2 


2.15 


X 


10 8 


1.04 


X 


10" 


-7 


4.24 


X 


10" 


-2 


3.56x10" 


-4 


1.31x10" 


-3 


236 tt 
92 u 


45.242 


12.404 


2si/ 2 


1.06 


X 


10 8 


1.76 


X 


10" 


-8 


8.47 


X 


10" 


-3 


1.60x10" 


-3 


2.00x10" 


-3 


236 tt 
92 u 


45.242 


12.698 


2^1/2 


3.02 


X 


10 9 


4.01 


X 


10" 


-7 


1.02 


X 


10" 


-2 


-1.26x10" 


-3 


1.27x10" 


-3 


236 tt 
92 u 


45.242 


16.871 


2p 3 /2 


3.10 


X 


10 9 


2.07 


X 


10" 


-7 


1.52 


X 


10" 


-2 


-9.86x10" 


-4 


5.01x10" 


-4 


238 tt 
92 u 


44.910 


12.073 


2^1/2 


1.11 


X 


10 8 


1.81 


X 


10" 


-8 


8.80 


X 


10" 


-3 


1.61x10" 


-3 


2.01x10" 


-3 


238 tt 
92 u 


44.910 


12.356 


2pi/ 2 


3.14 


X 


10 9 


4.17 


X 


10" 


-7 


1.06 


X 


10" 


-2 


-1.24x10" 


-3 


1.25x10" 


-3 


238 tt 
92 u 


44.910 


16.534 


2p 3 /2 


3.23 


X 


10 9 


2.16 


X 


10" 


-7 


1.56 


X 


10" 


-2 


-9.97x10" 


-4 


5.07x10" 


-4 


248 Cm 
96 Um 


43.380 


6.888 


2^1/2 


2.18 


X 


10 8 


3.25 


X 


10" 


-8 


1.78 


X 


10" 


-2 


1.92x10" 


-3 


2.16x10" 


-3 


248 Cm 
96 Um 


43.380 


7.190 


2^1/2 


5.47 


X 


10 9 


7.24 


X 


10" 


-7 


1.91 


X 


10" 


-2 


-5.96x10" 


-4 


5.99x10" 


-4 


248 frti 
96 Um 


43.380 


12.356 


2^3/2 


5.33 


X 


10 9 


3.54 


X 


10" 


-7 


2.20 


X 


10" 


-2 


-1.43x10" 


-3 


7.24x10" 


-4 



For the magnetic multipole transitions we consider the Ml transitions of the odd isotopes 
^ 5 Ho, ^ 3 Yb, |Mn, f 6 Fe, f° 9 K, IfGd, IfGd, If Re and ^ 7 Re. Numerical results for these ions 
are presented in Table|inj We present NEEC rates and resonance strengths with improved accuracy 
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TABLE II: Total RR cross sections for recombination into a given bound state of a bare ion, compared with 



results from Ref. 12011 . The nuclear excitation energy E cxc is given in the second column. The values from 
Ref. [20] are numerically interpolated by a spline routine to obtain the RR cross section at the resonance 
energy E c . 



isoiope 


-fiexclJiev; 


£/ c (Ke v j 


uruiiai 


i nis worK 


i?„f r9m 
ivei. [zuj 


66 u y 


I'X ^Q9 
1 j.jyL 


in 


is l/2 


8^9 


8^9 


68 tr 


78 <nQ1 


1 1 ^n 


i 'Sl/2 


7Q7 
ly 1 


7Q^ 

lyj 


174 V K 
70 Yb 


7£ AT 1 


A 8Q7 
4.07 / 


i -5l/2 


ZUoU 


9nen 

ZVoV 


154p j 
64 ua 


11'! (Y7 1 

iz j.v 1 1 


(sA nn^ 

04.UUJ 


J-Sl/2 


70 

ly 


70 

ly 


236 tt 
92 U 


45.242 


11 1 1 1 

11.113 


2S1/2 


245 


245 


236 tt 
92 U 


4J.Z4Z 


1 I.Ujo 


^Pl/2 


9Q^ 

zyj 


9Q4 
Zy4 


236 tt 
92 U 


4J.Z4Z 


1 J.OU1 


^>3/2 


99Q 

ZZy 


99Q 

ZZy 


238 tt 
92 U 


44 Q1 n 


1 n 789 
1U. / oZ 


^ s l/2 


9^9 

ZjZ 


9^3 


238 tt 
92 U 


44 Q1 n 


i n 7n^ 
iu. /uo 


4>l/2 




jUD 


238 tt 
92 U 


A A Q1 n 

44.7 iu 


1 C T£Q 

i j.zoy 


^3/2 




ZJU 


248 

96 Um 


4^ ^sn 


^ ^nn 


^ s l/2 


^4^ 
J4j 


^44 
J44 


IfCm 


43.380 


5.398 


2^1/2 


768 


769 


i 8 Cm 


43.380 


11.018 


2p 3 /2 


410 


410 


efHo 


94.700 


29.563 


1^1/2 


252 


252 


173 Yh 


78.647 


7.073 


l«l/2 


1410 


1410 


185 -Re 
75 ne 


125.358 


42.198 


!si/2 


212 


212 


75 ne 


134.243 


51.083 


lSl/2 


166 


166 


iMn 


125.949 


117.378 


lSl/2 


0.865 


0.849 


57 Fp 
26 re 


14.412 


5.135 


lSl/2 


216 


216 


40 K 
19^ 


29.829 


24.896 


1^1/2 


6.64 


6.55 



with respect to our previous results II 1411 . The electronic radial integrals are calculated numerically 
using the same type of wave functions for the bound and continuum electron as for the electric 
transitions. The reduced magnetic transition probability B (Ml) and the energies of the nuclear 
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levels are taken from Refs. [32, 



33, 



34j, 



35, 



36L131 



38, 



3H. 



Recombination into the K shell is 



possible for all the chosen ions, except for Gd. We also present results for recombination into the 
initially He-like ions of the gipGd and g| 7 Gd isotopes. 

TABLE III: Parameters of the NEEC total cross section and the interference term for various heavy ion 
collision systems involving magnetic dipole transitions. The notations are as defined in Table J] 



Isotope £ cxc (keV) ^(keV) Orbital 



Y n (l/s) T d (eV) S(beV) 1/Q f 



R u 



IfRo 

173 Y U 
70 

185 r>„ 
75 ne 

75 11(3 

|Mn 

57tp„ 
26 re 



155 



155 



Gd 
Gd 



IfGd 



94.700 
78.647 
125.358 
134.243 
125.949 
14.412 
29.829 
60.008 
60.008 
60.008 
54.533 
54.533 
54.533 



29.563 
7.073 
42.198 
51.083 
117.378 
5.135 
24.896 
45.784 
45.938 
46.722 
40.309 
40.463 
41.247 



Is 



1/2 



lSl/2 



Is 



Is 



Is 



Is 



1/2 
1/2 
1/2 
1/2 



2s x / 2 

2P1/2 
2^3/2 
2S1/2 
2P1/2 

2p 3 /2 



1.28 xlO 10 
7.32xl0 9 
2.62 xlO 10 
2.50xl0 10 
2.45 xlO 7 
6.21 xlO 6 
1.33xl0 7 
2.73 xlO 8 
2.40 xlO 7 
4.00 xlO 6 
4.16 x10 s 
3.68 xlO 7 
6.21xl0 6 



1.17xl0" 5 
4.80xl0" 6 
2.36xl0" 5 
2.47 xlO" 5 
1.75 xlO" 6 
2.56 xlO" 9 
9.47 xlO" 8 
1.97 xlO" 6 
1.86xl0" 6 
1.85xl0" 6 
4.37 xlO" 7 
2.71 xlO" 7 
2.56xl0" 7 



8.84xl0" 4 

1.26 

1.34 

1.16 
9.22 xlO" 4 
1.19xl0" 3 
2.27xl0" 3 
3.18xl0- 3 
2.94xl0" 4 
4.84xl0" 5 
2.86xl0" 2 
4.07 xlO" 3 
7.12xl0" 4 



-1.67xl0" 3 
-2.24 xlO" 3 
-2.58 xlO" 3 
-2.50xl0" 3 
-2.14xl0" 5 
-6.73 xlO" 5 
-1.46xl0" 5 
-1.25xl0" 4 
-1.85xl0" 5 
-1.81xl0" 5 
-1.25xl0" 4 
-2.00 xlO" 5 
-1.94xl0" 5 



2.90xl0" 3 
2.98 xlO" 3 
4.71 xlO" 3 
5.00xl0" 3 
3.10xl0" 3 
8.42 xlO" 5 
1.22 xlO" 4 
2.06 xlO" 3 
3.27 xlO" 3 
1.91xl0" 2 
3.00xl0" 4 
3.36xl0" 4 
1.82xl0" 3 



In Fig. [2] interference and scaled NEEC cross section terms are plotted as a function of the 
continuum electron energy for the Ml transition of 7 5 5 Re and E2 transition of 7 o 4 Yb, respec- 
tively. These are the isotopes with the largest values for the resonance strengths for the magnetic 
and electric multipole transitions, respectively. The NEEC cross section has the shape of a very 
narrow Lorentzian, with the width given by the natural width of the excited nuclear state, about 
2.4x 10" 5 eV for the case of 7 5 5 Re and 4.9x 10" 8 eV for the case of 7 Q 4 Yb. The interference term 
a int for both electric and magnetic cases is more than two orders of magnitude smaller than the 
NEEC terms cxneeo The magnitude of the interference term can be explained by investigating the 
contributions of the multipolarities that enter in the RR cross section <trr. While ctrr consists of 
an infinite sum of multipolarities, in the interference process only the RR photon with the multipo- 
larity of the nuclear transition participates. The main contribution to the RR cross section comes 
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FIG. 2: Interference and NEEC terms of the cross section for capture into bare -|| 5 Re (upper figure) and 
bare 7Q 4 Yb (lower figure) ions as a function of the continuum electron energy. The NEEC term is scaled by 
a factor of 10~ 2 . 



from the electric dipole El photon. The cross sections corresponding to the Ml and E2 photons 
are considerably smaller. In the case of yo 4 Yb, the E2 multipole accounts for only 121 b in the 
total RR cross section of 2080 b, while the Ml multipole for ^Re only contributes 0.5 b to the 
total RR cross section of 212 b. 

As an electron energy resolution in the order of 10 -5 eV and less can not be presently achieved 
in an experiment, we convolute the theoretical total cross section with the energy distribution of 
the electrons to give an orientation for measurements in near future. The energy distribution of the 
incoming electrons is assumed to be described by a Gaussian function with a width parameter s. 
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51100 



FIG. 3: The ratio R(E, s) in Eq. d67l ) for recombination into bare rhenium as a function of the energy of the 
continuum electron for three different experimental electron energy width parameters s. See text for further 
explanations. 

The RR cross section has a practically constant value on the energy interval of s. In order to 
demonstrate the magnitude of the NEEC and interference cross sections oneec an d er int compared 
to that of RR, we present in Fig. |3]the ratio of the convoluted cross sections, 



in the case of 75 Re as a function of the continuum electron energy for the three different exper- 
imental width parameters s = 0.5 eV, 1 eV and 10 eV. While for a width parameter s = 0.5 eV 
the contributions of the NEEC and interference terms can be clearly discerned from the RR back- 
ground, for presently more realistic widths in the order of eVs or tens of eV the values of the ratio 
R(E,s) are too small to be observed experimentally. 

IV. SUMMARY 

In this article we investigated the interference between NEEC and RR in an electron recom- 
bination process. We derived the interference cross section and expressed it with the help of the 
dimensionless Fano profile parameter. 

We calculated the interaction matrix elements for both electric and magnetic multipolarities 
using relativistic electronic wavefunctions. Nuclear excitations are described using a phenomeno- 



R{E,s) 



(Tneec(E,s) + a int (E,s 
a RR (E,s) 



(67) 
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logical nuclear collective model. The nuclear part of the matrix element is written by the help of 
the reduced nuclear transition probability whose value is taken from experimental works. For the 
quantization of the radiation field we use the multipole expansion. 

Numerical values for the Fano profile parameters and interference cross sections were obtained 
for various heavy-ion collision systems. The interference term in the total cross section of the 
recombination process is about two orders of magnitude smaller than the NEEC cross section. 
This is associated with the fact that from the infinite multipole expansion of the RR radiation, only 
the multipolarities corresponding to the type of nuclear transition interfere with the radiative decay 
photons following NEEC. The interference term has a narrow extent on the electron energy scale, 
which is related to the small natural width of the nuclear excited state. In order to simulate data of 
a recombination experiment, we convolute the total cross section with a Gaussian electron energy 
distribution of realistic width parameters. While for well-defined experimental electron energies 
the presence of NEEC could be discerned from the RR background, for larger width parameters 
both NEEC and the interference with RR become difficult to be observed experimentally. 

If the angular distribution of the emitted photons in the radiative decay of the nucleus following 
NEEC is different from the one of the RR photons, this can be used to identify the resonant 
process in the RR background. Calculations investigating a possible NEEC signature in the angular 
distribution of the emitted electrons are in progress. 
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